Thermoacoustic effects in supercritical fluids near the critical point: Resonance, 
piston effect, and acoustic emission and reflection 
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We present a general theory of thermoacoustic phenomena in supercritical fluids near the critical 
point in a one- dimensional cell. We take into account the effects of the heat conduction in the 
boundary walls and the bulk viscosity near the critical point. We introduce a coefficient Z{u) 
characterizing reflection of sound with frequency uj at the boundary. As applications, we examine 
the acoustic eigenmodes in the cell, the response to time-dependent perturbations, sound emission 
and reflection at the boundary. Resonance and rapid adiabatic changes are noteworthy. In these 
processes, the role of the thermal diffusion layers is enhanced near the critical point because of the 
strong critical divergence of the thermal expansion. 

PACS numbers: 05.70.Jk, 64.70.Fx, 62.60.-l-v, 65.40.De 



I. INTRODUCTION 

In highly compressible fluids, adiabatic changes take 
place with propagation of sounds and are much faster 
than the thermal diffusion When a one-component 
fluid is heated or cooled at a boundary, a thermal diffu- 
sion layer expands or shrinks to emit sounds, which then 
cause adiabatic changes in the interior (the thermal pis- 
ton effect). The density change in the boundary layer 
is enhanced near the gas-liquid critical point because of 
the strong critical growth of the isobaric thermal expan- 
sion. As a result, thermal equilibration times become 
shorter near the critical point at fixed volume, despite 
the fact that the thermal diffusion constant D tends to 
zero at the criticality 0, i, i, i, i, 0, i, i, 0, [H ■ K the 
boundary temperature is slightly changed, temperature 
homogenization occurs throughout the cell on a scale of 
the piston time 



(1.1) 



where L is the cell length, D is the thermal diffusion con- 
stant, and 7 = Cp/Cy is the specific-heat ratio growing 
near the critical point. See Appendix A for a derivation 
of Eq.(l.l). The time ti is much shorter than the isobaric 
equilibration time /AD by the factor (7 — 1)~^ <C 1. 
For example, we have ti = 6.3 x lO^e^'^^sec for CO2 
with L = 1cm. Hereafter, e = T/Tc — 1 is the reduced 
temperature near the critical point. The early experi- 
ments detected slow temperature and density changes in 
the interior on time scales much longer than the acoustic 
time L/c. The fast acoustic processes were examined by 
numerical simulations of the hydrodynamic equations of 
compressible fluids ITsj l . 

Ferrell and Hao found relevance of the heat con- 
duction in the boundary walls in transient heat transport. 
That is, the thermal boundary condition of a cell contain- 
ing a near-critical fluid crosses over from the isothermal 
to insulating one even for a metal boundary wall due to 
the critical divergence of the effusivity of the fluid p^ . 
The formula (1.1) should then be modified, because it 
is based on the isothermal boundary condition. More re- 



cently. Carles and Dadzie found that the bulk vis- 

cosity, which grows strongly near the critical point, can 
affect the hydrodynamics in the thermal diffusion layer. 
Gillis et al. [isj performed experiments of acoustic res- 
onance in xenon, where the frequency and attenuation 
of the resonating modes were measured. For such long 
wavelength sounds, the heat conduction at the boundary 
is the dominant damping mechanism relatively far from 
the critical point, while the viscous effect in the bulk be- 
comes more important closer to the critical point. They 
also presented thorough theoretical analysis of their data. 
The critical growth of the effusivity and the bulk vis- 
cosity of the fluid both serve to suppress the boundary 
damping, as confirmed experimentally and theoretically. 
Very recently, Miura et al. [l^, |20[ measured acous- 
tic density changes with precision of order 10~^g/cm'^ 
in near-critical CO2 on the acoustic time scale using a 
ultra-sensitive interferometer. They detected emission 
and traversal of sound pulses with width of order lO^sec, 
which were broadened as they moved through the cell 
and interacted with the boundary walls. Some of their 
data agreed with predictions, but most data remain unex- 
plained. Afterwords, part of the measured time-evolution 
of the density was numerically reproduced by Carles, ne- 
glecting the bulk viscosity [2]| . 

Some unique aspects of the supercritical hydrodynam- 
ics have been revealed by experiments [2^ . [23l . |23| and by 
simulations ^13, 25, 26, 27, 28, 29]. In Rayleigh-Benard 
convection, overall temperature changes are induced by 
plume arrivals at the boundary walls due to the piston 
effect, leading to overshoot behavior observed in experi- 
ments of '^He near its critical point [H, [1^, 113] ■ Signif- 
icant noises of the adiabatic temperature changes were 
predicted in turbulent convective states [1^, though not 
yet measured systematically. Recently three-dimensional 
simulations were performed |29j ]. Due to large thermal 
expansion in supercritical fluids, jet-like fluid flow has 
been observed around a heated boundary [2^ [2^. In 
these processes, the plume motions governed by the shear 
viscosity are strongly influenced by large thermal expan- 
sion around a heater and by rapid adiabatic density and 
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temperature changes achieved by sound propagation. 

In this paper, we aim to give detailed analysis of the 
linear hydrodynamics of supercritical fluids near the crit- 
ical point in a one-dimensional cell. We take into account 
the effects of the decreasing effusivity ratio and the 
growing bulk viscosity. In Section II, we will decompose 
fluid motions into sound modes and thermal diffusion 
modes with frequency lo. These two modes are mixed 
at the boundary under given boundary conditions, lead- 
ing to various thermoacoustic phenomena. In Section 
III, we will study the acoustic cigcnmodes determined to 
confirm the calculations by Gillis et a/.[l^ in the sim- 
pler one-dimensional geometry. We will also examine the 
response of the fluid to various time-dependent pertur- 
bations. Resonance is induced when the frequency of the 
perturbation is close to one of the eigenfrequencies, while 
nearly uniform adiabatic changes are caused in the inte- 
rior due to the piston effect at much lower frequencies. 
We will also examine sound emission and reflection at 
the boundary. In Appendix A, we will present a simple 
theory of the piston effect, which can be a starting point 
to understand the complicated calculations in the text. 
In Appendix B, the critical behavior of one-component 
fluids used in the text will be summarized. 



II. THEORETICAL BACKGROUND 



A. Linear hydrodynamics 



Here the prime denotes the differentiation with respect 
to X. We have two dissipative coefficients; one is the 
thermal conductivity A and the other is 



v, = {C + 4r?/3)/p, 



(2.4) 



where C, and 77 are the bulk and shear viscosities, respec- 
tively. Using the thermodynamic derivatives we may ex- 
press 5s and 5p in terms of 5p and 5T as 

pT5s = Cv[5T -b-^5pl (2.5) 
5p ^ -i-'^c^5p+{l~-f-^)as5T, (2.6) 

where — {dp/dp)s is the square of the sound velocity. 



7 = Cp/Cv 



(2.7) 



is the specific-heat ratio with Cp = pT{ds/dT)p and 
Cv = pT{ds/dT)p being the isobaric and constant- 
volume specific heat (per unit volume), respectively. To 
avoid cumbersome notation, we write 



dp\ , f \ -2 



(2.E 



For low-frequency sounds, the adiabatic relations Sp = 
GgST and Sp = bgST should hold. We use the following 
thermodynamic identities 



dp 
df 



] ={l-j-^)a,= pc^Cv/Tas 



(2.9) 




Near the critical point, we treat the hydrodynamic de- 
viations with spatial scales much exceeding the thermal 
correlation length ^, but the typical frequency lu can be 
higher than the relaxation rate of the critical fluctuations 
[soj . For such high frequencies ui, the bulk 
behaves as l/w, while ( (x e~^'^^ for ujt^ < 1 
, See Appendix B for more details. The other 
transport coefficients may be treated to be independent 
ofcj 1 . The critical singularity of the shear viscosity rj is 
negligible small, while the critical growth of the thermal 
conductivity A arises from the convective motions of the 
critical fluctuations taking place on a short time scale of 
order p^^/r?. Here we assume w <C v/p^^- 

The mass density, the temperature, the entropy (per 
unit mass), and the pressure are written as p, T, s, and 
p, respectively, with their small deviations being 6p, 5T, 
Ss, and Sp. The velocity in the x direction is written as 
V. These deviations depend on time t as exp{iLut) and 
vary in space along the x axis. We may assume lu > 
without loss of generality (see Eq.(3.1)). These devia- 
tions may be regarded as the Fourier transformations of 
the space-time dependent deviations with respect to time 
(= / die~"^*(- • •)). They obey the linear equations [s^ ]. 



iLuSp 
iLopv 
iojpTSs 



-pv , 

-Sp' -f pviv" , 
XST". 



(2.1) 
(2.2) 
(2.3) 



Next we consider small hydrodynamic deviations by 
assuming the space-dependence in the sinusoidal form 
exp(igx). From Eqs.(2.3) and (2.5) Sp and ST are related 

by 



ST 



(iuj + "/Dq'^)a, 



-Sp, 



(2.10) 



where D = X/Cp is the thermal diffusion constant. Equa- 
tions Eq.(2.1)-(2.3) give the dispersion equation between 
q and uj, 

= iwc2(l -7~i)q2. (2.11) 



If we set q = uj/cyX or A = (uj/cq)^, the dimensionless 
quantity X obeys the quadratic equation, 

X^-{l + A^+ jAt)X + (1 + 7A„) At = 0. (2.12) 

where we introduce two dimensionless coefficients repre- 
senting the dissipation strength [isj . 



At = itjjDjf?. 



(2.13) 
(2.14) 



If cij > 0, A^ and At are purely imaginary. The ratio 
A^/At — vijD grow strongly near the critical point (see 
Appendix B). 
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For given lu, Eq.(2.11) or Eq.(2.12) yields four solu- 
tions q = ±q± — ±uj/c^y X±, where X+ and X_ arc the 
solutions of Eq.(2.12) written as [l8!| 



X± = 1(1 



-7AtT^), 



2(1 + 7A^)At 
1 + A„+7At±S' 



where we define 



S = [(1 + A, - 7At)' + 4(7 - 1)At] 



1/2 



(2.15) 



(2.16) 



with ReS > 0. The second line of Eq.(2.15) follows from 
X+X- = (1 + 7A„)At. The modes with q — ±q^ repre- 
sent the sound, while those with q — ±g+ the thermal dif- 
fusion. We may define and (?+ such that Re{q^/u}) > 
and hnq^ < hold. It is convenient to introduce k and 
K by 



k = q^ = 



K = ig+ = 



(2.17) 



The argument of X- is in the range [0, 7r/2] for w > 0, 
leading to Imfc < 0, which implies that sound waves prop- 
agating in the positive x direction (cx e"*'^^) are damped 
with increasing x. 

As w ^ 0, we may treat A^, and A^ as small quantities. 
To their first order we find X+ ^ At and X_ = 1 -I- A^ -|- 
(7— 1)At so that K = y^iuj/D and k = u / c — iT sLo"^ /2c^ , 
where [Sji 



(C + 4,7/3)/p+(7-l)i? 



(2.18) 



is the attenuation constant in the long wavelength limit. 
We have |k| :s> \k\ at low frequencies. For example, 



10^ cm" 



1 and |fc| 
1 



10 cm ^ for a; = 10 s 



D = 10 ^ cm^s ^, and c = 10^ cm/s. In a cell with 
length L, the strength of the bulk dissipation of sounds 
is represented by the damping factor exp(— (5bL) with 



-Imfc 
,2 



r,cj72pc^ 



(2.19) 



where the second line is the low-frequency expression. 
Mathematically, we may consider the high frequency 
limit |A„| ^ 1 and |At| ^ 1 neglecting the frequency- 
dependence of the transport coefficients to derive the 
limiting behavior k — > [iiolvi)^!'^ and k — > {iu /^DY^'^ , 
though this limit is unrealistic. In this paper, we will 
assume |At| <C 7~^ in Eq.(2.34), because it is satisfied 
in realistic experimental conditions, as will be discussed. 



Solutions in a finite cell 



We consider small hydrodynamic perturbations behav- 
ing as e"' in a fiuid in a finite cell with length L. The 



density deviation can be expressed in the following linear 
combination, 



6p- 



ae 



Akx 



-j3e 



-ikx 



(2.20) 



The coefficients a, 6, a, and /3 depend on time as e'"*. 
The first and second terms represent the deviations in 
the thermal diffusion layers. The thickness of the layers 
is given by 1/|k|, which is assumed to be much shorter 
than the cell length L, so 



> l/L. 



(2.21) 



Then the second (first) term is virtually zero near x = Q 
[x = L). The third term in Eq.(2.20) represents a sound 
propagating in the negative x direction , while the fourth 
term a sound propagating in the positive x direction. 
From Eq.(2.1) the velocity is expressed as 



ILO , 



w = _iae - 6e'=(^-^)] - ^[ae**^^ - f3e-''''=]. (2.22) 
pK pk 

If the boundary walls are fixed in time, we should require 
w = at x = and L to obtain 



{K/ik){a — /3), 
-{K/ik)iae'''^ 



(2.23) 
(2.24) 



Note that the mass change in the thermal diffusion layers 
is {a + b)/K and that in the interior is a{e^''^ — l)/ik + 



/3(1- 



)/ik per unit area. From Eqs.(2.23) and (2.24) 



these two changes cancel, ensuring the mass conservation. 
Use of Eq.(2.10) gives the temperature deviation in the 
following linear combination. 



6T 



Vs 



H : 



ikx 



-ikx 



iu! — 'fDn^ 



-/Dk^ 



(2.25) 



Let Qo and Ql be the heat flux —XdT' at a: = and L, 
respectively. Use of Eqs.(2.23)-(2.25) gives 



/3 



P — ae 



2ikL 



bsh-1) 



rjikL 



Qi 



X{l + jAy) k^ + 



(2.26) 



where 63(7 - 1)/A = p/Ta^D with the aid of Eq.(2.9). 
The above quantities tend to the constant p/TagC in the 
low frequency limit. We may use Eq.(2.26) when Qo is a 
control parameter or when is measurable. 

The coefficients a, 6, a, and /3 can be determined if 
we specify the boundary conditions at a; = and L. 
Hereafter we assume no temperature discontinuity at the 
boundaries. In most theoretical calculations the bound- 
ary temperatures are fixed, but in some papers the bot- 
tom heat flux Qo is fixed [26j. In this paper, we consider 
a more realistic boundary condition of the temperature 
accounting for the thermal conduction in the boundary 
wall regions [ij, [11] . Here we assume that 5T tends to 
zero in the solid far from the boundaries without heat 
input. In the solid region (a; < 0), the temperature devi- 
ation then decays as (5T(0)e''"'^ with 



(2.27) 
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where A^, and Cw are the thermal conductivity and the 
heat capacity (per unit volume) of the solid, respectively. 
The 1/1 I is the thickness of the thermal diffusion layer 
in the solid and is assumed to be shorter than the thick- 
ness of the wall. Without temperature discontinuity at 
the boundary, the energy balance at a; = yields 



6T' — \yjK^ST/\ 

= a^{iuj/Df/'^5T, 



(2.28) 



where 5T and 5T' are the values at a; = 0. In the second 
line, the coefficient is the effusivity ratio p^. [isj. 



1/2 



(2.29) 



For CO2 in a Cu cell [2fl] we have = 3 x lO^e"-^^ The 
boundary temperature at a: = is fixed or 6T{0) — for 
Ou, — > 00, while the boundary is thermally insulating or 
[dST / dx)x=o = as Cu, — *■ 0. On the other hand, if the 
other boundary wall in the region x > L is made of the 
same material, the boundary condition at x = L reads 



ST' 



-a^{iuj/D)^^^ST, 



(2.30) 



with the same as in Eq.(2.28), where ST and 6T' are 
the values dX x = L. 

The boundary conditions at x = give Eqs.(2.23) and 
(2.28), from which we may readily calculate the reflec- 
tion factor Z = j3/a between the outgoing and incoming 
sound waves. It is convenient to introduce the combina- 
tion. 



W ■ 



a- (3 _ 1- Z 
a + (3 ~ 1 + Z' 



(2.31) 



because is a small quantity in our system. Some cal- 
culations yield a general expression. 



W = 



iui + 'fDJ'^ + \/iijjD{K'^ + k'^)/ayjK 



(2.32) 



in terms of k and k. In the case of a thermally insulating 
boundary, we have W — and Z — 1 by setting a^j — > 
in Eq.(2.32). The interaction of sounds and the boundary 
wall is characterized by Z or VF, where the wall proper- 
ties appear only through the effusivity ratio and the 
system length L does not appear. 



C. Adiabatic condition in the interior 

We will clarify an upper bound of the frequency, below 
which the sound motions in the interior are adiabatic or 
without entropy deviations. Under this adiabatic condi- 
tion, the results from the linear hydrodynamic equations 
can be much simplified. 

Far from the boundary walls or outside the thermal 
diffusion layers, we may neglect the localized modes 
to obtain the interior hydrodynamic deviations. From 



Eqs.(2.5), (2.6), and (2.25), those of the density, temper- 
ature, and pressure are related by 



Sp = [l+jDk^/iuj]bsST, 
Sp = [1 + Dk^/iLu]asST. 



(2.33) 



Here x and L — x are much longer than 1/|k|. The second 
terms in the brackets arise from a small entropy deviation 
in the interior. Since 7 > 1, the usual adiabatic relations 
hold in the interior imder the condition, 

\-fDk^/iuj\^j\AT\-^l or uj<s:c^/jD, (2.34) 

where At is defined by Eq.(2.14). This condition is well 
satisfied in the usual hydrodynamic processes. Even near 
the critical point, the time tad = jD/c^ remains very 
short. For example, tad = 7.6 x lQ-^'^e-°-^^ 

Under Eq.(2.34) we have S = 1 + A 
X_ in Eq.(2.15) are approximated as 



sec for CO2. 
so that and 



1 + 7A^ 
1 + A„ 



A. 



X^^l + Ay. 



The wave numbers k and k are expressed as 



1/2 



1 + A^ 
1 + 7A^ 



1/2 



k = 



uj/c 



VT+Ay 



(2.35) 



(2.36) 



We retain A^, , since it becomes appreciable near the crit- 
ical point because of the strong critical divergence of C,. 
As will be discussed in Appendix B, C — pc^Rsti, for 
Lut^ < 1, where Rb — 0.03 is a universal number and 
— / D is the characteristic time of the critical fluctu- 
ations with ^ being the correlation length. Carles found 
the dependence of k on the singular combination 7^ as in 
Eq.(2.36) [H, [i3|- By setting 7A„ = iujtB, we introduce 
a new characteristic time ts as 



(2.37) 



Then ts > once RbI > 1- For CO2, <b = 1-9 x 



10 



-15^-3.0 



sec. See Table 1 and Fig. 1 for the character- 



istic times with L = 1cm, where ts exceeds the acoustic 
time L/cior € < 3 x 10~^ and the modified piston time t[ 
(to be introduced in Eq.(3.16)) for e < 3 x 10^^. There 
can be a sizable frequency range with t^^ < uj < t^^ at 
small e, where k becomes independent of lu as 

K = {pc^hDCfl^ = {Rbi)-'^^C'- (2.38) 

The thickness of the thermal diffusion layer 1/\k\ remains 
longer than ^ by {Rb^)^^^- Also from the expression of k 
in Eq.(2.36) we write the sound dispersion relation as fc = 
llj/c*{uj), where we define the complex sound velocity [ij. 



'iu;)=c^l + A, 



(2.39) 



whose critical behavior will be discussed in Appendix B. 

As e is decreased, we first encounter the regime where 
W grows but ^ 1 and wis <C 1 still hold. However, 
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TABLE I: Parameters of CO2 in a Cu cell with L = 1 cm 
for e = 10~^ (first line), 10"" (second line), and 10"^ (third 
line). Times are in sec. 
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tf: X 10'' 


tB X 10^ 


L/c X 10* 


t\ 


t2 X 10® 


260 


5.0 


0.24 


1.9 


0.71 


1.0 


0.12 


3600 


0.63 


18 


1900 


0.83 


0.1 


1.7 


5 X 10* 


0.075 


1300 


1.7 X 10'' 


0.98 


0.08 


3.0 



10' 



10"' 
^10-* 

03 
0) 

1 10-s 



10-' 



10-' 









L7c 










t 2 


'^'"^"--^ 



10' 



10 " 10" 

e = T/T, - 1 



10'' 



FIG. 1: Characteristic times tf. in Eq.(B3), ts in Eq.(2.37), 
t2 in Eq.(2.43), A in Eq.(3.16), and L/c vs e = T /Tc - 1 for 
CO2 in a Cu cell with L = 1cm. 



the critical growth of W is eventually suppressed by the 
growing and C,. If we use Eqs.(2.35) and (2.36) under 
Eq.(2.34), we approximate W in Eq.(2.32) as 



W = 



(7- 



(1 + A„)X„ 



(2.40) 



where we define 



X„ = v/l+7A„+a-Vl + A„. (2.41) 

The limiting behaviors of are as follows: = 1 + 
for uj < fg" and X^, = (iutB)^^'^ for > t;g^(l + a~^)^. 
For w ^ tl^^, it follows the classical expression valid 



far from the critical point, 



(2.42) 



which is the result without the viscous effect and under 
the isothermal boundary condition. We may introduce a 
characteristic time ^2 defined by 



t2 = K,/(l + a.„)]2(7~l)2i?/c2 



(2.43) 



which includes the effect of the heat conduction in the 
wall. As shown in Table 1 and Fig. 1, t2 is very short 
even c omp ared with t^. In the literature (see Section 77 
of Ref.[3|), it is argued that the amplitude of a plane 
wave sound is decreased by the factor (7 — l)\/2Duj/c 
upon reflection at an isothermal boundary wall. This 
factor is obviously equal to 1 — |Z| = 2ReM^ if use is 
made of Eq.(2.42). 






' ' "1 ' — 1 ' — ^ 

E = 10"" 


1 ' — 












Re W 


10"' 










/\ 
\ 


10-" 


' , ' , 


' . 1_ 


10"^ 10"' 


1 



COL/JIC 



FIG. 2: ReW, ImW, and ]A„| = wC/pc^ vstjL/Trcat e = 10"^ 
(left) and 10"* (right) for CO2 in a Cu cell with L = 1cm. As 
LO exceeds tg^, KeW tends to saturate and ImW decreases, 
due to the growing bulk viscosity and the decreasing effusivity 
ratio. 



In Fig. 2, ReVl^ and ImW^ are displayed as functions of 
w at e = 10~^ and 10"''. While cutB < 1, they increase 
with increasing oj obeying Eq.(2.42). After oj exceeds 
vTTtS^ becomes y/Y^ in Eq.(2.41); then, ReM^ 
tends to saturate and ImM^ decreases. In fact, for cu ^ 
t^^(l + a~^)^, we have Xy = VtS^ and 
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(2.44) 



where Wq = 2.1 for CO2. Also as a function of e, 'ReW 
exhibits a maximum around the reduced temperature at 
which Ib ^ t^"^, as will be shown in Fig. 3. Growing 
a^^(cx e^^'^**) further serves to decrease W . Thus, even 
close to the critical point, we find |VF| ^ 1 and 



Z = 1 - 2T4^ 



(2.45) 



D. 



Hydrodynamic variables in the adiabatic 
condition 



Under Eq.(2.34) we obtain simple expressions of the 
deviations of the temperature, the pressure, and the en- 
tropy including A„. From Eqs.(2.5), (2.6), (2.9), and 
(2.20) we find 

ST - {l^^K)[5p\b + bs-\5p\,n (2.46) 



5s = [j-^^ (1 + A.)[M6, 
5p = -c^A.u[(5p]h + c^[(5/9]™: 



(2.47) 
(2.48) 



where [5p\b = ae"'^^ + 5e'*^^~^' is the density deviation 
localized near the boundaries and [Sp]in = ae^''^ + (3e~^'^^ 
is the interior density deviation. In deriving Eq.(2.46) use 
has been made of the thermodynamic relation (1 — 7)6^ = 
{dp/dT)p. Under Eq.(2.34), the entropy deviation Ss is 
localized near the boundaries, while the localized part 
of the pressure deviation Sp is nonvanishing to satisfy 
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iupv = —Sp' + pviv" = —5p' — iujviSp' — s- as x — > and 
L in Eq.(2.2). 

We examine the deviations close to the boundary at 
X = by asssuming Eqs.(2.23) and (2.28) and setting 
Za ~ (3. In this case the density ratio [5p\b/[5p]in tends 
to (7 — + 7^i> as X 0, so that 



bsST = (a + /3) 



1 



I - —x/l + jA^e- 



1 



(7~1)A 



, (2.49) 
, (2.50) 



where < x <^ 1/1^1- Note that the second terms 
in the brackets in Eqs.(2.49) and (2.50) tend to unity 
for ujtB ^ (1 + a~^)~^, which can be the case of very 
large C. In the original work the pressure homo- 
geneity and the isobaric relations among the hydrody- 
namic variables were assumed in the thermal diffusion 
layers. We recognize that the pressure homogeneity and 
the isobaric condition hold only in the low frequency limit 
a;«tBi(l + a-i)-2. 



III. APPLICATIONS 
A. Acoustic modes in a cell 

Gillis et al. [l^ calculated the acoustic eigenmodes 
for their experimental geometry, taking into account the 
growing and ^. In the following, we will present a 
simpler version in a one-dimensional cell, < x < L, 
taking into account these two ingredients. In this case 
Lu is treated as one of the eigenvalues and is complex, 
while we have assumed w > in the previous section. 
Then uj should have a positive imaginary part for the 
stability of the system. Here Viw = (1 -f i)y/ujj2 for 
UJ > 0, while ^/iu; = (1 - i)^/\uJ\/2 for w < 0. The latter 
follows from the requirement that the real part of k in 
Eq.(2.17) should be positive. For general complex lu, all 
the quantities introduced so far should be functions of 
u! analytic for Re(ia;) > or for Imuj < 0. Therefore, 
X±(-u;) =X±(w*)* and 



Z(-cj*) = Z(w)*, W{-uj*) = W{ujy 



(3.1) 



where is the complex conjugate of lo. 

Under the boundary conditions (2.28) and (2.30), the 
interior density deviation is expressed as 

Sp = aie''"'-" + Ze-'^'-'), 

= a'ie''''-^-''^ + Ze-'''^^-'=^), (3.2) 



in terms oi Z = (3 /a. The first and second lines follow 
from Eqs.(2.28) and (2.30), respectively, and should co- 
incide so that a'e^^^ — (3 and a'Ze^'*^^ — a, leading to 
2-g-*fci — Z~^e^^^. We now find the condition of the 
eigenmodes. 
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n = 1 ■ 

























1 ' 

1 n = 4 ■ 























E = T/T^ - 1 



10"^ 10 10'' 

E=T/T^-1 



FIG. 3: Normalized damping constants a\, a^, and a\ -\- 
defined by (3.5) and (3.6) vs e = T /Tc - 1 for n = 1 (left) 
and n = A (right). The resonant frequencies are close to 
ui = mxc/L ~ n X 4 X lO^'sec"^, which are exceeded by t~j^^ 
close to the critical point as marked by the arrows. 



where -I- corresponds to even modes and — to odd modes. 
Namely, the density and temperature deviations are even 
(odd) functions of x — L/2 for the the even (odd) modes. 
For W = {1 - Z)/{1 + Z) calculated in Eq.(2.32) or 
Eq.(2.40), we obtain 



W ^ -itim{kL/2) 
= icot{kL/2) 



(even modes) 

(odd modes). (3.4) 



Since W is small, the eigenfrequencies cun [n = 
1, 2, 3 • • •) are nearly equal to mrc*/L where c* is defined 
by Eq.(2.39). In the case |IF| ^ 1 the leading correction 
from W can be written as 

w„ = {mT + 2iW~\ )c*/L, 

= [1 + iax + ia(^ + ■ ■ ■]Reu;n, (3.5) 

where n = 2, 4, • • • for the even modes and n = 1, 3, • • • for 
the odd modes. We assume small bulk damping | Ai,| <C 1 
and define 

ax = —ReW, = hA,\ ^ ^^i^T^ (3-6) 
rnr 2 zpcL 

where ax represents the boundary damping and the 
bulk damping. The resonance quality factor [l^ 
is equal to 2 (a a + 0;^) in our notation. The resonance 
frequency including the shift is given by the real part. 



Reujn = (riTT - 2ImT4^ + • • ■)Rec*/L. 



(3.7) 



±e 



ikL 



(3.3) 



The frequency u in c* and W may be equated with 
Reujn = n-KcjL. 

In Fig. 3, we show ax, a^, and the sum ax + as 
functions of e in the regime cjtj < 1 for the odd mode 
of n = 1 and the even mode of n = 4 for CO2 in a 
Cu cell with L — 1cm. We notice the following, (i) 
For such long wavelength sounds, the boundary damp- 
ing is relevant far from the critical point, but the bulk 
damping eventually dominates close to the critical point, 
(ii) In accord with the discussion around Eq.(2.44), ax 
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decreases on approaching the criticaHty in the region of 
ts > uj~^ ^ L/mrc, with a maximum at Lotg ~ 1. As a 
result, the curve of the sum ax + is flattened consid- 
erably around Ib ~ uj~^. 

These theoretical results are consistent with the ex- 
perimental data by Gillis et a/. [3]. They performed the 
resonance experiment over a wide range of tot^ (up to 
about 200) to measure the frequency-dependent bulk vis- 
cosity. In agreement with the theory 0,0], or uC,/ pc^ 
became independent of e in the high-frequency regime 
ijjt^ > 1 (see Eq.(B5) in Appendix B). 



B. Periodic perturbations 

Periodic perturbations may be applied to a fluid in 
a cell in various manners. Resonance can occur when 
the frequency lj is close to Rew„. It is sharp for small 
IroLOn- We will give three boundary conditions at a: = 
leading to resonance. We assume the boundary condition 
Eq.(2.30) at X = L. Then use of Eq.(3.2) yields a = 
f3Ze~^^^^ . The interior density deviation is of the form, 

5p = pe-'^'-" + fiZe''"''-^'''^ , (3.8) 

where the term proportional to Z arises from the reflec- 
tion at X = L. 

The bulk damping of the reflected waves is represented 
by |e~*'^^| — er^"^ . From Eq.(2.19) 5b is expressed as 

6b= AB{uj/^cf, (3.9) 

in the low frequency regime < 1. We find Ab — 0.5 x 
lO^'^cm and 0.03cm at e = 10^'^ and 10^**, respectively, 
for CO2. In the relatively high frequency range uj > 
(2As)-i/27rc/L, the factor e-^'''^{(x e'^^"^) becomes 
negligibly small. Then, near the boundary, dp consists 
of the outgoing wave only, resulting in no resonance. On 
the other hand, in the high frequency regime Lot^ > 1, 
Eq.(B6) gives 

dB=(^lmAy/2c^0.27uj/TTc, (3.10) 

which means that a sound emitted at x = reaches the 
other end with the damping factor e~° -^^ for the first 
resonance frequency lo = ttc/L. 



1. Temperature oscillation 

In the first example, the temperature in the wall region 
X < is oscillated, while the boundary walls are mechan- 
ically fixed. More precisely, we require ST ^ oc e*'^* 
as X ^ -00; then, ST{x) = e''™^[5r(0) - T^] + in the 
region x < 0. The thermal boundary condition at x = 
is then given by 

ST' = a^iiuj/D)^/\6T-T^), (3.11) 




10 I I I I I L. 

1 2 3 4 5 

COL/TtC 

10"^ r 1 1 1 1 r 




10 '■^ I 1 1 1 1 1 

1 2 3 4 5 

COL/JtC 

FIG. 4: Absolute value of the response function -Rt(i^) in 
Eq.(3.12) vs ljL/ttc on a semi-logarithmic scale for e = 10"'^ 
(upper panel) and 10~^ (lower panel) applicable for CO2 in a 
Cu cell with L — 1 cm. 

as a generalization of Eq.(2.28). Some calculations using 
Eqs.(2.23) and (2.25) give the response function defined 
by Rt = l3/bsTw in the form, 

Rt = lil+JDky^cu)^^0^. (3.12) 

Notice that Rt diverges as Rt = cW /2iL{u! ~ ujn) for 
UJ = LUn (n = 1,2,---) in the complex lo plane from 
Eq.(3.3). Under the adiabatic condition Eq.(2.34), the 
interior temperature deviation is expressed as 

ST = b-^Sp = (e-'*^'^' + Ze'^''-'^'^^)RTT^. (3.13) 

Furthermore, neglecting jDk'^ /iuj in Eq.(3.12) and using 
\W\ < 1 (see Fig. 2), we obtain 

Rt'^W/[1- {l-4W)e-^'''^]. (3.14) 

In Fig. 4, we plot the absolute value \Rt\ calculated 
from Eq.(3.12) vs the normalized frequency ldL/'kc at 
e = 10~^ and 10"'', using the data for CO2 in a Cu cell 
with L = 1cm [20]. It exhibits peaks at w = mrc/L as 
expected, but its peak heights do not exceeds 1/2 due to 
the small factor 1 — Z = 2W in the numerator in Eq.(3.12) 
As discussed below Eq.(3.9), the resonant peaks should 
disappear for ujL/ttc > (2v4b)"^/^, where we may neglect 
g-2ifcL obtain Rt = W. These results are in 

accord with Fig. 4, since {2Ab)^^^^ =30 and 4 for e = 
10~^ and 10~*, respectively. 

In the low frequency case cu <^ c/ L, the interior devi- 
ations become nearly homogeneous. Figure 1 indicates 
that ts can much exceed L/c very close to the critical 
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point, while |A„| <C 1 holds. Thus, retaining 7A„, we 
set e*''^ ^ 1 + ifcL and 1 + A„ 5^ 1 and use Eqs.(2.40) 
and Eq.(3.12) to find 



1, 



(3.15) 



where tx is defined by Eq.(l.l) and by Eq.(2.41). If 
<^ tg^, we further have Rt = {\/iiot\ + l)~^/4 with 



A = (1 



% = (1 



(3.16) 



which is related to in Eq.(2.43) by t[t2 = L'^/Ac'^. The 
t'l first decreases as ti ^ e^'^^ for a^, 1 but finally 
weakly increases as tia^ = g-0.22 ^^j. ^ -[_. See Fig. 1 
for the curve of t[ . As will be discussed in Subsection III 
C, t[ is the piston time including the effect of the wall 
heat conduction [l3 |. 

We need to know when \y/iujtiXy\ ^ 1 holds. It holds 
for CO 1/t'i under the condition, 

(1 + a-'yhB/t[ = (1 + a-')-HB/ti « 1. (3.17) 

If a„ < 1 for CO2 in a Cu cell, the above condition 
becomes a^tB/h = 2 x IQ-^e'^ '^'^ / L"^ < 1 with L in cm 
, which is well satisfied for e ^ 10^^ with L = 1cm. If 
LL)t[ < 1 under Eq.(3.17), we find 

ST - (3.18) 

in the interior. Note that the reverse condition of 
Eq.(3.17), d^tB/ti > 1, holds extremely close to the crit- 
ical point, where \^/iujtiXy\ > 1 and Rt — l/Aiuj\/titB 
are obtained for uj ^ (ii^s)"^^^. See the discussion be- 
low Eq.(3.32) for the relaxation behavior in this ultimate 
regime. 

Zhong et al. [Tl| measured a density change induced 
by boundary temperature oscillation in near-critical "^He, 
where the frequency was very low (w/27r < 2Hz) and the 
bulk viscosity was not important. However, they could 
measure in-phase and out-of-phase response in agreement 
with the original theory [3]. 



2. Mechanical oscillation 

In the second example, the boundary wall at x = is 
mechanically oscillated without heat input from outside. 
This is the case in the usual acoustic experiments using 
a piezoelectric transducer [3]. Let Uu,(cx e*"*) be the 
applied displacement amplitude; then. 



at X in Eq.(2.12) 
Eq.(3.8) we obtain 



(3.19) 

Assuming Eq.(2.30) and using 



P = - 



kRt 



1 — ^Dk^/ilu 
(7 - DXyRTpu^ 



where the first line is general and the second line is the 
approximation under the adiabatic condition Eq.(2.34). 
Since the response is proportional to Rt, resonance oc- 
curs as in the previous case of temperature oscillation. 

In the low frequency case lu <C c/L, the interior density 
change is nearly homogeneous and 



1 



1 



^/iujtiXy + 1 



L 



(3.21) 



which is the counterpart of Eq.(3.14). As discussed be- 
low Eq.Eq.(3.16), ^/uJtiXy is large in the relatively high- 
frequency range uj ^ ^/t'l under Eq.(3.17). Thus the 
interior density deviation behaves as 

Sp ^ pu^,/L {l/t[ < w < c/L) 

^ ViuJhXypu^,/L (w<l/t;), (3.22) 

under Eq.(3.17). The volume change mostly occurs in 
the bulk region for 1/t'i ^ u; ^ c/ L and in the thermal 
diffusion layers for uit'i <^ 1. 



3. Heat flux oscillation 



10 











e = 10" 


3 
















FIG. 5: Absolute value of the response function Rq{uj) in 
Eq.(3.24) vs LdL/nc on a semi-logarithmic scale for e — 10~^ 
(upper panel) and 10~^ (lower panel) applicable for CO2 in a 
Cu cell with L = 1 cm. 

In the third example, we apply a heat flux Qq = 
—X{dT/dx)x=o oc e'"* at a; = assuming the boundary 
condition (2.30). It is convenient to introduce a dimen- 
sionless response function Rq by 



(3.20) 



(3.23) 
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Then Eqs.(2.26) and Eq.(3.8) give 



Rq 



ick/D 



(1 + 7A^)(A:2 + k2)(1 - Ze-2''fci) 
1 

(l + A^)3/2(l-Ze-2*fei)' 



(3.24) 



where the first hne is general and the second line holds 
under the adiabatic condition (2.34). In the complex 
uj plane, Rq has poles w^, which are equal to uj2n in 
Eq.(3.5) with system length changed to 2L. Thus Rq 
grows for w = mrc/L for w < {2Ab)~^^^ttc/ L. In Fig. 5, 
we plot the absolute value \Rq\ as a function of uL/ttc 
for e = lO^-^ and 10"'*. We can see that |i?Q| is larger 
than \Rt\ in Fig. 4 roughly by two orders of magnitude. 

The behavior of Rq in the low frequency range co <C 
c/L is very different from that of Rt, however. From the 
second line of Eq.(3.24) we have Rq = l/{2ikL + 2W) to 
obtain the counterpart of Eq.(3.15), 



Ro = 



2^JTuA2 



[1 + 2^fiZA[X^\ 



(3.25) 



Under Eq.(3.17) we find that Rq = (l + a^*)/2V«tI^ for 
LU < l/t[ and Rq = l/AiLUy^t^ for l/t[ < w < c/L. 

In this situation we may calculate the heat flux Ql at 
X = L. From Eq.(2.26) it is written as 



Qi 



(1 - Z)e- 



-ikL 



1- Ze 



-2ikL 



(3.26) 



which vanishes for Z = \ (or for — 0) and becomes 
small with increasing 5bL. Near the resonance frequency 
m:c/L, the ratio Ql/Qo behaves as {—l)^W/[i{ujL/c — 
mi) + mTC,/2pc^ + W\. The low frequency behavior for 
Lo <^ c/L is given by 



Ql = (1 + 2V^X.„)"iQo. 



(3.27) 



From the discussion below Eq.(3.16), we find Ql — Qo 
for Lut'i <C 1 under Eq.(3.17). That is, an applied heat 
flux passes through a near-critical fluid on the time scale 
of t[ under Eq.(3.17), due to the piston effect. 



C. Thermal and mechanical piston effects 

L Boundary temperature change 

In the original papers of the piston effect [3*], the 
boundary temperatures at x = and L were both raised 
by a common small amount Ti at t = 0. Subsequently, 
the boundary temperatures were held flxed for t > 0. In 
this paper, we examine the effects of finite [3, 
and large C, [H, We suppose that the system was 
in equilibrium for t < and the temperature in the 
wall region x < was instantaneously raised by Ti at 
t = without external heat input in the other wall re- 
gion X > L. The boundary conditions are then given by 




FIG. 6; Relaxation function *(*) in Eq.(3.29) defined by 
Eq.(3.30) vs t/t{ at e = 10'^ 10"", and 10"^ for CO2 in 
a Cu cell with L = 1 cm, which is applicable for t ^ L/c. 
The time t'l is defined in Eq.(3.16). Functional form of ^'(f) 
as a function of t/t'i is rather insensitive to e since Eq.(3.17) 
holds for this case. 



6T{x, i) — !■ Ti as a; — > —00 and 6T{x, i) — s- as x~L 00 
for t > 0. All the deviations vanish for t < 0. 

The Fourier transformation of the interior temperature 
deviation 6T{x,t) with respect to t is given by Eq.(3.13) 



with T„ = Txe'^^/iuj (since dte- 
inverse Fourier transformation gives 



5T{x,t) 



27ria; 



-ikx 



Ze 



ikx — 2ikL 



l/iuj). The 



)Rt, (3.28) 



where the integration is in the range [— (xi,oo]. Under 
Eq.(2.34), W = W{lo) and Rt = i?r(w) are given by 
Eqs.(2.40) and (3.12), respectively. The integrand is 
analytic (without singularities) in the lower half plane 
Imw < and hence the integral is nonvanishing only for 
t > 0. 

In the time region t 3> L/c we may neglect the space 
dependence of 5T{x, t) in the interior and use the simple 
expression (3.15) for Rrit). It then follows 



6T{t) = Ti*(t)/2, 



(3.29) 



where we introduce the dimensionless relaxation function 
"^{t). Its Fourier transformation reads 



1 



iuj{\/TujtiXv + 1) 



(3.30) 



The inverse Fourier transformation of the right hand side 
of Eq.(3.30) may be transformed into an integral along 
the positive imaginary axis Ivauj > 0. With being 
defined by Eq.(2.41), we generally find ^'(i) > for t > 0, 

^'(i) ^ t/y/ttt^ as t ^ 0, and ^{t) = 1 - {t'Jnty^^ H 

as i — > 00. In particular, not very close to the critical 
point, we may neglect the bulk viscosity and take the 
limit t_B ^ 0; then, X^ ^ + and *(t) ^0(3), 
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where ^'o(s) is a universal function of s = t/t[ expressed 
as [| 

r°° du e""'' 

*o(s) = l-/ ^f— -l-e^erfc(yi), (3.31) 

where erfc = 1 — ere is the complementary error function 
and = 2(s/7r)i/2 for s < 1 and *o = 1 - 2(7rs)~i/2 
for s 3> 1. 

In Fig. 6, we display '^{t) as a function of t/t'i at 
e = 10-3, 10-^ and 10"^ for CO2 in a Cu cell with 
L = 1 cm. For e = 10"^ we can see *(t) = \E'o(t/ti), 
where is/t'i ~ 0.02 from Table 1. The discussion below 
Eq.(3.16) indicates that \l/(<) approaches unity on the 
time scale of t[ as long as Eq.(3.17) is satisfied. This 
is the case even for e — 10^^, where ts/t'i = 21 from 
Table 1. In fact, if is/i'i ^ 1 and ^ 1, we may 
set ^/HJtiXy = iuj^titB + a~^^/liJt^, where the second 
term is relevant in '^{t) under Eq.(3.17), again leading 
to '^{t) = 5'o(t/ti) for t > a^ts- However, the reverse 
condition of Eq.(3.17) holds extremely close to the critical 
point, where Rt = l/[iuj^/titB + 1] holds yielding [l3| 

^!{t) = 1 - cxp(-t/ViTtB). (3.32) 

The new relaxation time \/titB here grows as \/titB ^ 
1.0 X 10~''Le~°-^^sec for near-critical CO2. 

Assuming the isothermal boundary — 00), Carles 
and Dadzie examined the bulk viscosity effect in the 
thermal equilibration [l7l |. Their relaxation function is 
obtained if we set Xy = {1 + iojtBY^^ in Eq.(3.30). 
Then a new viscous regime appears for ii 3> is with 
'i'{t) being given by Eq.(3.32), while the usual piston 
regime is encountered for ti ^ tB- For CO2 we have 

tB/h ^ 2.7 X 10-20L-2g-4.63^ ^ j^^l^g 

e = 0.6 X lO^'' with L — 1cm. In our calculations 
based on Eq.(3.17), the different predictions have arisen 
from the reduced temperature dependence of or the 
crossover of the boundary condition into the insulating 
one. 



2. Volume change 

We suppose a volume change by moving the boundary 
wall at X = by a small length ui instantaneously at 
i = pj]. We assume the thermal boundary conditions 
(2.28) and (2.30) at x = and L. As in Eq.(3.26), the 
complete interior density deviation is the inverse Fourier 
transformation of Eq.(3.8), where /? is given by Eq.(3.20) 
with = uie'^'^*- 1 iuj . 

Here we are interested in the late stage Ljc^ where 
the interior deviations depend only on t. The inverse 
Fourier transformation of Eq.(3.21) gives the interior de- 
viations, 

hp{t) = h,6T{t) = [1 - ^(t)\pu^lL, (3.33) 

where ^(t) defined by Eq.(3.30) represents the effect of 
the thermal diffusion layers at a; = and L. The above 



form with = was derived in Ref.[l|. If u\ > 0, 
the interior is adiabatically heated by b~^pui/L on the 
acoustic time scale L/c after the volume change, while the 
boundary wall temperature is almost unchanged. Subse- 
quently, the thermal diffusion layers become effective as 
reverse pistons and the interior temperature deviation 
decays as (1^11)^1'^ . 

The reverse piston effect itself generally occurs on the 
time scale of t\ after a near-critical fluid was adiabatically 
heated or cooled. Miura et al. observed such a process 
after a pulse-like heat input (see Fig. 2 in Ref.|20l|). 

D. Emission of sound 

We examine sound emission at the boundary at a; = 0. 
We neglect the incoming wave reflected at the other end 
X — L and consider the semi-infinite limit L 00. 

The problem is simple in the case of boundary wall 
motion. An emitted sound propagates with the velocity 
c and integration of the continuity equation gives the 
density deviation, 

5p{x,t) = pvi(t — X / c) / c, (3.34) 

where vi it) is the velocity of the boundary. The localized 
part of the density deviation (the term proportional to 
a in Eq.(2.20)) should be small when differentiated with 
respect to time. In fact, under the adiabatic condition 
(2.34), Eqs.(3.12) and (3.20) lead to 

^ - 2(l^+f„)c ^"^"- ^^-^^^ 

for 6^2*'=^ -> 0. If we set H- A„ ^ 1 and Z ^ 1, 
the above relation becomes (3 = iupuwlc, leading to 
Eq.(3.34). Here iuuw is the Fourier transformation of 
vi{t) multiplied by e''^*. Thus Eq.(3.34) holds on time 
scales longer than (even when the time scale of vi{t) 
is shorter than Ib)- 

A sound is also emitted when a time-dependent heat 
flux Qo{t) is supplied at the boundary at x = 0. From 
Eq.(3.8) the Fourier transformation of the interior den- 
sity deviation Sp{x, t) is of the form /3e''^^ with /3 be- 
ing given by Eq.(3.23). Under the adiabatic condition 
Eq.(2.34) we may use the second line of Eq.(3.24) to find 
the convolution relation, 

Sp{x,t) = ^(^^^ y* dT^x,t^T)Qo{T). (3.36) 

The memory function <I>(x,t) is defined for t > as 

$(x,t)= / ^e'-*-''=-(l + A,)-3/2, (3.37) 
J 27r 

where A^, — iujC./pc^ — iujRBt^ with Rb = 0.03. 
The time integration of this function is normalized as 
dt^{x,t) = 1. From the integration in the region 
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u) < t^^ we obtain the long-time behavior $(0,t) = 
(4</7ri^)i/2e"*/*': with = Rst^ at a; = 0. Since this re- 
laxation is rapid, we may set <&(0,<) = 6{t) ((5-function) 
at a; = on time scales longer than or when QQ(t) 
varies slower than t^. Furthermore, if the distance x is 
not large such that the bulk damping is negligible in the 
region < a; < L, we may set ^{x, t) = S{t — x/c) to find 
the simple formula for the emitted sound, 

6pix, t) = ^ — j Qo{t - x/c), (3.38) 

as the counterpart of Eq.(3.34). On the other hand, use 
of Eq.(B6) for A„ gives the short-time behavior. 



(3.39) 



valid in the time region t with d/2i> = 0.088 (see 
Appendix B). This behavior is detectable only for an in- 
crease of Qo within a time shorter than . 

Miura et al. applied a stepwise heat flux with Qo = 
0.183 X 10^ to find a stepwise outgoing sound with Sp/ p = 
2.2 X 10^^ for CO2, where Qo is in cgs units (erg/cm^sec) 
[20j . Our theoretical expression (3.38) becomes Sp/p = 
1.38 X IQ-^^Qo with the aid of {dT/dp)s = T^/QMpa 
for CO2 fssj . For their experimental Qo our theory gives 
5p/ p = 2.55 X 10^^ in fair agreement with the observed 
density change. Furthermore, they could generate sound 
pulses with duration of order 10/isec by applying short- 
time heat input. They were interested in the adiabati- 
cally increased energy i?ad = P J dx5p{x, t)/ p in the pulse 
region per unit area. Here Eq.(3.38) yields 0] 



^ad-f (^ ) Q, 



(3.40) 



where Q = / dtQQ{t) is the total heat supplied. The ratio 
Es^d/Q represents the efficiency of transforming applied 
heat to mechanical work. Theoretically, it is given by 
{dT/dp)sp/T as in Eq.(3.40land is equal to 1/6.98 = 
0.14 for near-critical CO2 [35|. The measured values of 
the ratio E^d/Q were in the range 0.11 — 0.12 again in 
fair agreement with our theory. 



E. Reflection of sound 

Reflection of plane wave sounds is discussed for an 
isothermal boundary in the textbook of Landau-Lifshitz 
[33 |. Miura et al. [20] observed reflected pulses passing 
through a detector in the cell. Their shapes gradually 
flattened after many traversals within the cell, resulting 
in the interior temperature homogenization. At present, 
it is not clear how to understand their data. Here, as a 
first step, we will derive some fundamental relations on 
sound reflection. 

We consider a pulse approaching to the boundary at 
a; = in the semi-infinite limit L 00. Reflection takes 



place upon its encounter with the wall. The density devi- 
ations of the incoming and outgoing pulses are obtained 
as the inverse Fourier transformation of Eq.(2.20). Ne- 
glecting the bulk damping in the neighborhood of the 
boundary, we may express them as pi{t + x / c) and Po{t' 



x/c), respectively. Using a(w) 
(3 = Za, we obtain 



/ dre-"^'' p^{T) and 



Po{t) 



dhJ 

2^ 



dt'Z{uj)e"^^'-'"^p^{t') 



(3.41) 



The interior density deviation is the sum 5p{x, t) — pi{t-\- 
x/c) + Po{t — x/c). Since Z(0) = 1, the excess mass is 
invariant upon reflection as 



dtp,{t) = / dtpo{t) 



(3.42) 



This relation holds if we integrate a long tail of the re- 
flected pulse Po{t) at large t (see Eq.(3.47)). 

If Pi{t) changes much slower than t^, we may set 
Z'^l-2W with W = (7- l)\/iA^/X„ from Eq.(2.40). 
In this approximation we may rewrite Eq.(3.40) in the 
following convolution form. 



Po{t) = p.(<)- / dTx{T)[p^{t-T)- p,{t)] 

Jo 



P^{t) 



dTx{T)Pt{t-T) 



(3.43) 



where x{t) — dx{t)/dt and pi{t) — dpi{t)/dt. From 
Eq.(3.41) the function xiP) is the inverse Fourier trans- 
formation of (1 — Z)/iui = 2W/iuj. Some calculations (in 
the complex lj plane) give x{^) in the integral form. 



xit) 



dn 



Re 



-Sis 



(3.44) 



where s — t/ts is the scaled time, Re[- • •] denotes taking 
the real part, and Vl — ^ = — 1 for > 1. The 
dimensionless parameter Sr is defined by 



2(7-1)/d7^ 



(3.45) 



which decreases near the critical point as £r = 4.3e°'^^ 
for CO2. The function xit) depends only on s and a^,. 
For ^ 1 we have x(0/^>- — e^''/^/o(s/2) with Iq 
being the modified Bessel function, while for a^i ^ 1 
we have x{t)/£r = 1 — $o(s/a^) with $0 being defined 
in Eq.(3.31). Thus xi^) changes on the scale of t'g = 
tB{l + a^^)^^ and its limiting behaviors are as follows: 



= (l + a-i)-i(^,s)-i/2 + ... (,s^oo), 

^ l-2a-i(s/7r)i/2 + ... (s->0). (3.46) 



In addition, the second term of Eq. (3.43) representing the 
distortion is negative (positive) when Pi{t) is increasing 
(decreasing). This initial drop is because of heating and 
expansion of the pulse at the boundary. 
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Prom the first line we obtain x{t) — — (i2/7r)^/^t~'^/^ 
for t'g. (i) Let pi{t) is peaked in the region |t| < to; 
then, for t ^ t'g and to, the first line of Eq.(3.43) gives a 
long-time tail of the reflected wave, 



[p(i)]t,ail - AM(t2/7r)l/2i-3/2 



(3.47) 



where AM is defined by Eq.(3.42). If to > t'^, the total 
mass behind the peak is given by the time integral of 
the tail Eq.(3.47) in the region [to,oo]. Thus the mass 
fraction behind the peak is (4t2/7i'io)^''^- For CO2 this 
quantity is estimated as 10~^e~°'^^t(5^^^^ with Iq in sec for 
ttw ^ 1. (ii) As another example, we consider a stepwise 
change, where pi{t) is equal to for t < and to a 
constant pi for t > with to being the transient time. 
Then, for t 3> t'^ and to, the second line of Eq.(3.43) 
gives a longer tail, 



[p(t)]tail=Pl(4t2A)'/'i-'/^ 



(3.48) 



The bulk viscosity does not appear in these tails. 

When pi{t) changes much slower than t'^, only the long 
time behavior of x{t) is relevant in Po{t)- From Eq.(3.43) 
we find the following convolution relations. 



Po{t) = P.{t) + \ - 




dT 
7372 

■ dT 



[p,{t - t) - p,{t)] 



p^{t-T), (3.49) 



from which the long-time tails (3.47) and (3.48) read- 
ily follow. The above expressions contain only ^2 in 
Eq.(2.43) and not ts- They are widely applicable far 
from the critical point (where ts becomes short). With 
decreasing e for the isothermal boundary (a^j > 1), t2 
grows and the distortion of the reflected pulse increases as 
long as the pulse width is longer than ts- However, if the 
pulse width is shorter than t'g, the distortion decreases 
on approaching the critical point since e,- in Eq.(3.45) 
decreases. 

As a simple illustration, let us consider a Gaussian 
pulse pi{t) = piexp{—t'^/2to), where pi is the pulse 
height and to is the pulse width. Since its Fourier trans- 
formation is (27r)^/^pitoe~'^ *a/^^ we may readily calcu- 
late Po(t). In Fig. 7, we plot the normalized pulse defor- 
mation defined by 



Fait) = [poit) - p^{t)]/ipiy/h/to)- 



(3.50) 



The curve (a) is for the limiting case ts/to — > and 
flu, = 00, while ts/to = 10 and a^, — 0.63 in (b), and 
ts/ta = 50 and a^, = 0.63 in (c). In Table 1 we have 
— 0.63 and ts — 1.9msec at e = 10^** for CO2 in a 
Cu cell, where pulses with to ^ts are well possible [20| . 
We recognize that the distortion is negative for t <to and 
is positive for t >to (in accord with the comment below 
Eq.(3.43)) and that the distortion is decreased as ts/to 
is increased or for shorter pulses due to the bulk viscosity 
growth. 




t/lr 



FIG. 7: Scaled pulse deformation Fait) defined by Eq.(3.50) 
vs scaled time t/to for a Gaussian incoming pulse with width 
to for ts/to — » and am = 00 (a), ts/to = 10 and a™ — 0.63 
(b), and ts /to = 50 and a™ = 0.63 (c). 



IV. SUMMARY AND REMARKS 

In summary, we have examined various thermoa- 
coustic effects in one-component supercritical fluids in 
a one-dimensional geometry. We summarize our main 
results. 

(i) In the linear hydrodynamics, sound modes and ther- 
mal diffusion modes arc both present as in Eq.(2.20), 
depending on given boundary conditions. The latter 
modes can be absent only for the insulating boundary 
condition = 0. The calculations are straightforward 
and the final expressions are much simplified under 
the adiabatic condition (2.34) or for low frequencies 
uj ^ c^/jD. It is remarkable that the bulk viscosity 

appears in the combination uj^C,/ pc^ — ujtB as first 
pointed out by Carles [l^, [l3| • The resultant character- 
istic time ts grows as e~'^ °, while the life time of the 
critical fluctuations t^ grows as e~^^^. 

(ii) We have introduced the reflection factor Z{uj) as the 
ratio between outgoing and incoming sounds. Using Z 
or = {1 — Z) / {1 + Z) we have examined the acoustic 
eigenmodes, the response of the fluid to applied oscilla- 
tion of the boundary temperature, the boundary heat 
flux, and the boundary position. To these thermal and 
mechanical perturbations, resonance is induced when 
the frequency of the perturbation is close to one of the 
eigenfrequencies, while nearly uniform adiabatic changes 
are caused in the interior at much lower frequencies 
owing to the piston effect. 

(iii) We have also examned the response to a stepwise 
change of the boundary temperature and the boundary 
position. The relaxation time is given by the modified 
piston time t[ in Eq.(3.16) first introduced by Ferrell 
and Hao [l^. It is equal to the original piston time ti 
in Eq.(l.l) for the isothermal boundary ^ 1 and to 
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a~'^ti for the insulating boundary a„j ^ 1. 

(iv) As the critical point is approached, the role of the 
thermal diffusion layers is eventually diminished both 
by decreasing of the effusivity ratio a^, and by growing 
of the bulk viscosity while the bulk sound attenuation 
becomes increasing stronger. The bulk viscosity effect 
in the thermal diffusion layer is thus masked by its 
enhanced effect in the bulk. 

(v) For CO2 in a Cu cell, the boundary becomes 
thermally insulating for e ^ 10^^. This suppresses the 
bulk viscosity effect in the thermal diffusion layers as 
long as Eq.(3.17) holds or for e > 10^^. In this case, the 
viscous regime predicted by Carles and Dadzie emerges 
for e < 10~^ [3, To increase this crossover reduced 
temperature, the cell length L needs to be shorter. For 
the wall materials in Ref. |18|. this crossover occurs much 
closer to the critical point p^ . 

(vi) We have also examined sound emission and reflection 
at the boundary, which are elementary hydrodynamic 
processes but seem to have not been well examined 
|34| . For emission, the formulas (3.34) and (3.38) 
are valid for a mechanical piston and a thermal heat 
input on time scales longer than t^. For reflection, 
Eq.(3.43) with Eq.(3.44) holds on time scales longer 
than t^. The formula (3.49) is the classical one valid on 
time scales much longer than ts, where the distortion 
of the outgoing pulse increases on approaching the 
critical point. For pulses shorter than tB, the distortion 
of the outgoing pulse is decreased as can be seen in Fig. 7. 

In this paper, we have treated near-critical fluids in 
one phase states. However, more challenging are hydro- 
dynamic effects in in two phase states, where latent heat 
transport, wetting dynamics, and Marangoni convection 
come into play in addition to the piston effect [36l[37l. [38t . 
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Appendix A: Simple theory of the piston effect 



the interior and are given by 



5p- 



6p p (dT\ SQ 



T\dp 



(Al) 



where L — V/A is the cell length. The interior tem- 
perature deviation is caused adiabatically as ST = 
{dT/dp)sSp and is written as 



= (7 - l)6Q/CpL, 



(A2) 



where Cp = pT{ds/dT)p is the isobaric specific heat and 
use is made of Eq.(2.9). If the boundary temperature is 
raised by Ti at t = 0, we have 5Q ~ Cp£{t)Ti in the 
early stage, where £(t) = y/Dt is the thickness of the 
thermal diffusion layer. If we set ST = Ti/2, Eq.(l.l) is 
reproduced. 

The relation (Al) also follows from our formula 
Eq.(3.38). Let the heat input rate Q{t) from the bound- 
ary to the fluid change slowly compared to the acoustic 
time ta = L/c. We suppose a time interval with width 
St^ta, in which Q{t) is almost unchanged. Since St/ta 
is the traversal number much larger than unity, the adia- 
batic pressure and density increases in the interior region 
are given by 



Sp. 



5p St p fdT 



ta cT \ dp 



(A3) 



as a result of superposition of many steps. In terms 
of the incremental heat supply SQ = QSt we reproduce 
Eq.(Al). 

Appendix B: Summary of critical behavior 

Let a one-component fluid be on the critical isochore 
{p — Pc) with small positive e = T /Tc — 1 near the gas- 
liquid critical point. The physical parameters used in 
the table 1 and the figures are given below. Hereafter 
v{= 0.63), 7(= 1.24), and a{= 0.10) are the usual critical 
exponents. Data of near-critical CO2 can be found in 
Refs.fiO,!!!. 

Our hydrodynamic description is valid when the spa- 
tial scale under investigation is longer than the correla- 
tion length i = ^oe"^, where ^0 = l-SA for CO2. The 
constant- volume specific heat Cv = pT{ds/dT)p and the 
isobaric specific heat Cp = pT{ds/dT)p are expressed as 



Av[e''^ + B], Cp 



(Bl) 



Here we give a simple derivation of ii in Eq.(l.l). Let 
us apply a small heat SQ to a fluid from the bound- 
ary per unit area in the one dimensional geometry. 
The volume expansion of the thermal diffusion layers is 
given by {dT/dp)sASQ{t)/T, where A is the area of the 
heater surface and use is made of the Maxwell relation 
{dp~^ /ds)p = {dT/dp)s- The interior density change 
and the pressure change Sp are nearly homogeneous in 



For CO2 on the critical isochore, the coefficients are given 
by Av = 26.3kBn*, B = 0.9, and Ap = 2.58kBn*, where 
n* — pc/kBTc = 1.76 X lO^^cm"^. The specific-heat ratio 
7 grows strongly as 7oe^'''^" if the background (cx B) is 
neglected, where 70 = 0.1 for CO2. The sound velocity 
and the constant- volume specific heat are weakly singular 
as cx e" (if the background is neglected l|). We have 
set c = 2.3 X lO^eO "^ cm sec"! for CO2 
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The thermal conductivity A grows such that the ther- 
mal diffusion constant D behaves as 

D = \/Cp = kBT/GmjC = Dae^ , (B2) 

where Dq = 4.0x IQ-W^sec-^ for CO2. Thus A oc e^""*". 
The relaxation time of the critical fluctuations with size 
^ increases as 

k = ^VD = toe-^', (B3) 

where to = 0.56 x lO^^^sec for CO2. The shear viscosity 
77 is only weakly singular and may be treated as a con- 
stant independent of e and lo to make rough estimates. 
However, the zero-frequency bulk viscosity C grows very 
strongly as 

C^Pc'RbM, (B4) 

where Rr is a universal number estimated to be about 
0.03 [ll,[3|. For CO2, C/p = 0.9 x lO-^e-^+s^cm^sec-i, 
so C/pD = Ay/ At = 0.02e-2-«>+2a (ggg Eqs.(3.13) and 
(3.14)). In the high frequency regime ujt^ ^ 1, the com- 
plex sound velocity in Eg. (2. 39) becomes asymptotically 
independent of e [3ll. [33. l33j . Thus, 

c*{lo) ^ c{iojt^f'^r (B5) 



Since the exponent al&v is small, we may set A^, = 
iLoC^I pc? = (ci/3i>) In(iwi^). Thus, in this high frequency 
regime, ImA„ tends to the following universal number, 

ImA„ = 7rd/6z> ^ 0.27 x 2/7r (B6) 

In the high frequency regime ujt^ > 1, Ay remains to 
be as a small quantity and the frequency-dependent bulk 
viscosity defined by ({cu) = pc^A^/iuj decays roughly as 
l/iw with increasing uj. 

Furthermore, in our thermoacoustic problems, we have 
introduced the time tB in Eq.(2.37), which behaves as 

tB = 4e-^^-^+", (B7) 

where = 1.7 x 10"^^sec for C02. The effusivity ratio 
in Eq.(2.29) decreases as 

= aO e^--^/^. (B8) 

For a^, ^ 1 the boundary wall crosses over from an 
isothermal one to an thermally insulating one on ap- 
proaching the critical point. For example, between Cu 
and CO2, we have aJJ, = 3 x 10'^ where a^, < 1 is 
reached for e < 1.6 x 10""*. The was smaller for the 
walls used in Ref.[il[li]. 
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